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Abstract 



We evaluate regions of parameter space in the minimal supergravity model 
where "unbounded from below" (UFB) or charge or color breaking minima 
(CCB) occur. Our analysis includes the most important terms from the 1- 
loop effective potential. We note a peculiar discontinuity of results depending 
on how renormalization group improvement is performed: One case leads 
to a UFB potential throughout the model parameter space, while the other 
typically agrees quite well with similar calculations performed using only the 
tree level potential. We compare our results with constraints from cosmology 
and naturalness and find a preferred region of parameter space which implies 
m-g < 725 GeV, rrig < 650 GeV, < 225 GeV and m^^ < 220 GeV. We 

discuss the consequences of our results for supersymmetry searches at various 
colliding beam facilities. 
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I. INTRODUCTION 



The Minimal Supersymmetric Standard Model (MSSM) is one of the leading candidate 
models [|l| for physics beyond the Standard Model (SM). In this theory, one begins with 
the SM particles (but with two Higgs doublets to ultimately ensure anomaly cancellation); 
supersymmetrization then leads to partner particles for each SM particle which differ by spin- 
|. Supersymmetry breaking is implemented by adding explicit soft supersymmetry breaking 
terms. This procedure leads to a particle physics model with ^ 100 free parameters, which 
ought to be valid at the weak scale. 

To reduce the number of free parameters, one needs a theory of how the soft SUSY 
breaking terms arise, i.e., how supersymmetry is broken. In the minimal supergravity model 
(SUGRA) [0, supersymmetry is spontaneously broken via a hidden sector field vacuum 
expectation value (VEV), and the SUSY breaking is communicated to the visible sector via 
gravitational interactions. For a flat Kahler metric and common gauge kinetic functions, this 
leads to a common scalar mass mo, a common gaugino mass mi/2, and common trilinear 
and bilinear terms Aq and Bq at some high scale Mqut — Mpi, where the former choice 
is usually taken due to apparent gauge coupling unification at ~ 2 x 10^^ GeV. The high 
scale mass terms and couplings are then linked to weak-scale values via renormalization 
group evolution. Electroweak symmetry breaking [Q], which is hidden at high scales, is 
then induced by the large top-quark Yukawa coupling, which drives one of the Higgs field 
masses to a negative value. Minimization of the scalar potential allows one to effectively 
replace B by tan (3 and express the magnitude (but not the sign) of the Higgsino mass /i in 
terms of M^. The resulting parameter space of this model is thus usually given by the set 
{mo,mi/2, Ao,taxi P, and sign(/i)). 

Not all values of the above 4 + 1 dimensional parameter space of the minimal SUGRA 
model are allowed. For instance, the top and bottom Yukawa couplings are driven to infinity 
somewhere between Mz and Mgut for tan/3 < 1.5 and ^ 50 (depending on the value of rrit) 
For other parameter choices, the lightest chargino or lightest slepton (or top squark) 
can be the lightest SUSY particle, which would violate limits on, for instance, heavy exotic 
nuclei. For yet other parameter choices, electroweak symmetry breaking leads to the wrong 
value of Mz, so these parameter choices are ruled out. In addition, there are cosmological 
bounds from the relic density of neutralinos produced in the Big Bang Requiring the 
universe to be older than ~ 10 billion years leads to only a subset of the parameter space 
being allowed, although this bound could be evaded by allowing for a small amount of R- 
parity violation. Finally, certain regions of parameter space are rejected by negative searches 
for sparticles at colliding beam experiments, such as those at LEP and the Fermilab Tevatron 

i- 

An additional constraint on minimal SUGRA parameters can be obtained by requiring 
that the global minimum of the scalar potential is indeed the minimum that leads to ap- 
propriate electroweak symmetry breaking. In the SM, there is only a single direction in 
the field space of the scalar potential, so appropriate electroweak symmetry breaking can 
be assured. For the MSSM, the plethora of new scalar fields which are introduced leads to 
many possible directions in field space where minima could develop which are deeper than 
the standard minimum. Thus, parameter choices which lead to deeper minima should be 
excluded as well, since they would lead to a universe with a non-standard vacuum. 
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Constraints along the preceeding lines were developed in Refs. in the early 1980's, 
using the renormalization group improved tree level effective potential. It was noted by 
Gamberini et. al. ^ that the renormalization group improved tree potential was subject to 
large variations due to uncertainty in the correct scale choice Q at which it was evaluated. 
Inclusion of 1-loop corrections served to ameliorate this condition. Recently, Casas, Lleyda, 
and Munoz have made a systematic survey of all possible dangerous directions in scalar 
field space that can potentially lead to minima deeper than the standard one. These have 
been categorized as field directions that are either unbounded from below (at tree level) 
(UFB) or that lead to charge or color-breaking (CCB) minima. For simplicity, their analysis 
uses the tree-level scalar potential but evaluated at an optimized mass scale where 1-loop 
corrections ought to be only a small effect. Working within the minimal SUGRA model, 
they considered models with Bq = Aq — itlq or Bq = 2mo and showed that significant regions 
of parameter space could be excluded via this method. 

In the present work, one of our main goals is to delineate the parameter space regions 
where non-standard potential minima develop in such a manner as to facilitate comparisons 
with other constraints, including recently calculated results on the neutralino relic density 



T0| and parameter space regions favored by fine-tuning considerations In addition 



expectations for supersymmetry at LEP2 |T2[, the Tevatron MI and TeV33 upgrades []T3 
the LHC and NLC [15| have been calculated within the minimal SUGRA framework. 
We also compare the non-standard vacuum constraints with the various collider expectations 
and draw some conclusions. For instance, combining the non-standard vacuum constraints 
with the most favored parameter space regions from fine-tuning and cosmology suggests that 
the Fermilab TeV33 upgrade stands a high chance to discover SUSY via the W1Z2 3£ 
signal! 

In the present work, we also adopt a somewhat different calculational scheme from that 
employed in Ref . . For all field directions considered, we implement renormalization group 
(RG) improvement to calculate the 1-loop effective potential. We find that the inclusion 
of the 1-loop correction has important consequences. The 1-loop correction almost always 
represents a significant contribution to the tree level potential. Nevertheless, our overall 
results agree very well with those of for a "proper" choice in RG improvement scheme. 
For reference, we shall call this the "a-case," and it represents our main results. However, 
for other choices of RG improvement, we find that the 1-loop correction can be so dominant 
as to lead to unbounded from below (UFB) potentials everywhere in parameter space: We 
refer to this as the "cj-case." Because this is a multi-scale problem, it is not entirely clear 
how to proceed with RG improvement, and it is this ambiguity that leads to the two cases 
above. The validity of our results hinge on two main assumptions. The first concerns the 
adequacy of cutting the expansion at 1-loop. It is beyond the present analysis to ascertain 
the significance of the 2-loop contribution in any of the cases considered. However, because 
of the dependence of our results on the details of RG improvement, we believe that the 
2-loop contribution may be important. Secondly, we have only included the contribution of 
the top-stop sector in our calculations of the 1-loop correction. It remains to be determined 
if the inclusion of the other fields will significantly affect the results. 

We note briefly that the work of |]T6[ and [|1^] advances the idea that we may indeed exist 
in a false vacuum and that the tunnelling rate from our present vacuum to a UFB or CCB 
vacuum might be small relative to the age of the universe. In this case, the following derived 
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constraints would not be meaningful. Such a philosophy must, however, be reconciled |T8 
with the fact that we live in a world in which the cosmological constant either vanishes or 
is extremely small. This is empirical, albeit indirect, evidence for some mechanism which 
seeks to enforce the principle that "the cosmological constant of the true vacuum is zero." It 
is difficult to conceive of circumstances where we could tenably entertain both the idea that 
we are living in a false vacuum and the idea that the smallness of the cosmological constant 
has a natural solution, since this would require a principle which would set the cosmological 
constant to zero in a false, broken vacuum while simultaneously leaving the true, broken 
vacuum with a large negative cosmological constant. 

The organization of the paper is as follows. In Section II, we review the MSSM scalar 
potential and give a brief summary of the UFB and CCB directions delineated in Ref. [§. 
In Section III, we present our calculational procedure and in Sec. IV present results of 
our scans over SUGRA parameter space. In Section V we give a brief summary of our 
results. Detailed formulae for the effective potential in various UFB and CCB directions 
are included in Appendix A, while some computationally useful formulae for evaluating the 
effective potential in the limit of large VEVs is presented in Appendix B. 



II. DANGEROUS DIRECTIONS IN FIELD SPACE 



The scalar potential of the MSSM can be written as 

V = V^ + Vj, + Vsoft , 

where 
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where r = (ri, T2, r^) are the SU{2) Pauh matrices, and A = (Ai, . . . , Ag) are the Gell-Mann 
SU{?>) matrices. 
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^soft = $Z^Jj0ar + {B/lK'^d + C.C.) 

a 

+ J2 {Au^Vu.Ui'^uQi + Ad^VdA^dQi + A^.y^.ei^dLi + c.c.) , (2.4) 

i 

and the supcrpotential W is 

= I] (l/«i«i*«<3i + VdA^dQi + VeA^dLi) + nK^d ■ (2.5) 

i 

In the above, 0^ runs over the scalar components of the chiral superfields, and a, i are gauge 
group and generation indices, respectively. Qi (Lj) are the scalar partners of the quark 
(lepton) SU{2)l doublets, and TZj, di, and are the scalar partners of the SU{2)l singlets. 
$u and are the two Higgs doublets. When all the above summations are performed, one 
is left with a very lengthy expression for the scalar potential. In the following, usually only 
a small number of scalar fields develop VEVs, so only a subset of the many terms of the 
scalar potential are relevant. 

For the usual breaking of electroweak symmetry in the MSSM, only and $d develop 
VEVs, so that the relevant part of the above potential is just 

Vo = m?|$rf|2 + m^|$„|2 + m^($„$rf + h.c.) 
where the masses appearing above are defined as 

2 2,2 

= m^^ + n , 

2 2,2 
^2 = ^<S>u + ' 

The 1-loop contribution to the scalar potential is given by 

where is the field dependent squared mass matrix of the model, and rrii is the mass of the 
i*'* particle of spin Sj. In the 1-loop correction, we shall always include only the contribution 
from the top-stop sector, this being generally the most significant. 

Minimization, at a scale Q usually taken to be Mz, yields two conditions on the param- 
eters 

1 2 ml — ml tan^ (3 
2^^ ~ tan2 13-1 

where m| = g'^v'^ — {g'"^ -\- gDv"^ /2, v"^ — vl + v% and 

Bii = ]^{ml + ml) sm2(5 , (2.12) 
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(2.6) 



(2.7) 
(2.8) 
(2.9) 



(2.11) 



with tan (3 = Vu/vd, and where the barred masses are the 1-loop analogs of ( |2.7| - |2.8|) . At this 
point, the minimum of the tree-potential is 



min 



- {m\ — m\) + {m\ + cos 2/5 



(2.13) 



The field dependent top and stop masses are given by 




li,2 = + 2 ^"^Q3 + "^i) + 4"^l 



± Vi2^'"^«~"^^3^^T2^^"^'^~^"'^^''°^^^J +"^? (A + /icot/3)2 , (2.15) 




where = g^v"^ /2. We will compare the value of the potential at the MSSM minimum 
with the value at the minimum for other field configurations and use this to reject MSSM 
scenarios with false vacua. 

The dangerous directions in field space have been categorized in Ref. |Q as various UFB 
and CCB directions. For the UFB directions, the trilinear scalar terms are unimportant. 
To find the deepest directions in field space, one searches for directions where the D-terms 
of Eq. ( |2.3| ) will be small or vanishing. The various UFB directions are characterized as 

• UFB-1: Here, only the fields and obtain VEVs, with < >=< $rf > in order 
to cancel D-terms in ^.6| . 

• UFB-2: In addition to the VEVs < > and < <l>d >, one has a VEV for the 3rd 
generation slepton field in the u direction: < L3 >y. The VEVs are related as in 
Eq. (A.3). 

• UFB-3a: In this case, the relevant VEVs are <$„>,< L3 >e=< 63 >^, and 

< L2 >y. This direction reputedly leads to the most stringent bounds on parameter 
space. The VEVs are related as in Eq. (A. 10). 

• UFB-3b: This case is similar to UFB-3a, but instead of the first two slepton fields, 

< Q3 >1=< ds >2 develop VEVs. The VEVs are related as in Eq. (A. 16). 

The various CCB directions each involve a particular trilinear coupling. For each trilinear 
coupling, there are two relevant directions: CCB(a) (equivalent to Casas et. al. CCB-1), 
and CCB(b) (which combines the CCB-2 and CCB-3 cases of Casas et. al.). The CCB(a) 
direction is not relevant for the top trilinear term. Summing over the various trilinear terms 
and CCB directions can yield at least 17 possible directions (some other possible directions 
lead to essentially the same constraints). For illustration, we investigated the following cases, 
which include cases with the largest and smallest Yukawa couplings. 

• CCB(a)-UP: The relevant VEVs are <$„>,< Qi >„, < ui >, < Q3 >1=< c/3 >^, 
and < L3 with < $d >= 0. The VEVs are related as in Eq. (A. 22). 

• CCB(b)-UP: The relevant VEVs are <$„>,< ^d >, < Qi >«, < Ui >, and 

< L3 >y. The VEVs are related as in Eq. (A.28). 
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• CCB(b)-TOP: The relevant VEVs are $d >, < Q3 >„, < M3 >, and 

< L3 The VEVs are related as in Eq. (A.42). 

• CCB(a)-ELECTRON: The relevant VEVs are < <^d >, < ^1 >e, < ei >, and 

< Q.^ >l=< Us >^ with < >= 0. The VEVs are related as in Eq. (A. 58). 

• CCB(b)-ELECTRON: Lastly, the relevant VEVs are < $rf >, < >, < Li >e, 
and < ci >. The VEVs are related as in Eq. (A. 63). 



III. CALCULATIONAL DETAILS 

The standard procedure for studying the MSSM has been, as summarized above, to fix 
the parameters Bq and fiQ in order to achieve symmetry breaking as dictated by ( pTf ) and 
( ^.81 ) with V = 174 GeV. Furthermore, the choice in minimization scale being in the Mz 
range is dictated by the desired vacuum expectation value and is validated by the use of 
the 1-loop correction. Unfortunately, this method does not lend itself to the present task, 
since a priori the minimum of the potential in a given configuration is unknown. The task 
is complicated since we must be able to probe the potential for significantly different field 
values. 

In order to validate the use of the 1-loop effective potential, one must ensure that not 
only the couplings be perturbative but that the logarithms be small as well. This is the 
process of RG improvement [jl9[. In problems with only one mass scale, RG improvement 



is straightforward. The logarithm appearing in the 1-loop correction can be made small, 
indeed to vanish, for any choice in field value by an appropriate choice in renormalization 
scale Q. This procedure yields the Q-independent, 1-loop RG improved potential. For the 
cases in which we are interested, there are several mass scales. Since in general no scale 
exists that simultaneously makes all the logarithms vanish, we settle for the scale at which 
the logarithms are simultaneously, optimally small. In this way, we construct the 1-loop 
effective potential. Note that in such cases the 1-loop correction does not vanish and indeed 
may represent a significant contribution to the tree level part. In our subsequent results, we 
always include the 1-loop correction in our evaluation of the effective potential. Figure 1 
demonstrates the significance of the 1-loop correction for a representative case with Aq = 0, 
mo = 100 GeV, mi/2 = 200 GeV, tan/?(M^) = 2, and < 0. In this example, we have 
employed the a-scheme (see below) for RG improvement. From Fig. 16 we see that the 
difference between the value of the tree potential at the minimum and the 1-loop effective 
potential in the UFB-3(a) direction is almost a factor of four. There is also an effect in the 
standard MSSM direction as seen in Fig. la. This particular point in parameter space is 
ruled out since Kn™"^'''' < Knin^^- implement the RG improvement procedure as follows: 

(1) At each RG scale Q find the field value that minimizes the function 

/(0,Q) = E[iogH(0,Q)/Q'}-x]' , (3.1) 

i 

(2) store this value and the corresponding Vi(0, Q). This results in the function (0, Vrgi{(P)) 
whose minimum can then be calculated. 

Given the standard numerical procedure involved in RG studies, in which Runge-Kutta 
routines are used to integrate over Q, the above procedure of finding the optimal (p at each 
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scale Q is the most efficient, since we construct the RG improved potential simultaneously 
as we evolve the RG parameters of the MSSM. Also, we note that including coefficients 
in Q) (as they appear in AVi) leads to pathological results. Namely, at all scales Q, the 
/-minimizing field value tends to zero. This is a pathology of the method we are employing. 
Had we instead fixed cf) and found the /-minimizing value of Q, this problem would not 
be present. Because of the ambiguities in RG improvement in multi-scale problems, we 
examined several prescriptions for constructing the effective potential. Eq. (2.9) with x = 
3/2 led to the a-case. We also tried x = 0. We used the functional form / = /n^{m^/Q^}, 
and tried the top and stop masses as possible choices for m. All of these choices led to the 
u;-case. In this case, the potentials along the UFB-3(b) and CCB(a)-UP directions were 
unbounded from below everywhere in parameter space. Although this result is interesting, 
we believe that 2-loop leading logarithms may remedy this curious situation. Therefore we 
discount these results for the moment. Along all other directions considered, the tu-results 
were similar to the a-results as one expects. Figure 2 displays results in the cj-case similar 
to Fig. 1. Comparing Figs. 1 and 2, the 1-loop effective potentials are essentially identical 
for the MSSM direction. And although the UFB-3(a) potentials are clearly different, the 
uj potential in this direction remains well-behaved. Figure 3 displays the logarithm of the 
VEV versus x two vacuum directions. In the UFB-3(a) direction, there is some minor 
change in the VEV but not very significant as x varies from the a to uj value. In constrast, 
the VEV suffers a drastic, discontinuous jump at around x = 1 in the CCB(a)-UP direction. 
This same discontinuity occurs in the UFB-3(b) case. 

From Fig. 3, we see that altering x (which effectively changes our Q choice) changes our 
results dramatically in the CCB(a)-UP and UFB-3(b) directions. These two directions are 
special in that their 1-loop contributions are dominant in the large VEV domain. This comes 
in particular from contributions to the stop masses by the yl^<P terms (see (A. 19-20,25-26)); 
bear in mind that the quadratic (or G±) contributions in x are always small or zero. For 
these cases, in the large VEV domain, the 1-loop correction, AVi, is obviously unstable 
against Q. For the case depicted in Fig. 3, we have verified that Vi(Q, = 10^^ GeV) has 
a large, negative slope and changes sign at Q ~ 4 x 10^^ GeV, leading to a potential which 
is very unstable against variations in scale choice. We assume this instability is due to a 
need to include higher terms in the effective potential for these cases. 

It is interesting to note that the standard procedure for computing the MSSM minimum 
gives results that differ, sometimes greatly, from the above RG improved procedure. We 
find that the 1-loop correction to the MSSM potential does a poor job of stabilizing the 
potential against Q near Mz- This was alluded to in The result is surprising given 

the conventional lore that the 1-loop correction stabilizes the potential in the electroweak 
range. A more detailed study of this issue is in progress. In Fig. 4, we demonstrate the 
problem by displaying the Q evolution of the VEV using two different methods in the case 
Aq = 0, mo = mi/2 = 100 GeV, ta.iaP{Mz) = 2, and fi < 0. The solid line represents the 
evolution as dictated by the RG gamma functions of the Higgs fields (see |2l|). The dashed 



line represents the tracking of the minimum of the 1-loop effective potential. Here we have 
included the contributions from all particles in AVi. 

For all cases, we have explored the dangerous directions delineated in Ref. 0; these 
directions were obtained using only the tree level potential. A better procedure would involve 
optimization of the full 1-loop effective potential. This may in principle lead to even more 
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precipitous directions. However, since this procedure leads to very unwieldy expressions, 
we have opted to explore the tree-level-derived directions in field space, although in these 
directions we make comparisons using the 1-loop corrected scalar potential. 

We emphasize the importance of using the 1-loop correction to the scalar potential, since 
its inclusion can alter the depth of the minimum significantly as was evident in Fig. lb. 
Furthermore, the startling results of the w-case were a consequence of using the 1-loop 
correction to compute the potential. Had the 1-loop correction been ignored, all choices 
of RG improvement we tried would have led to results similar to ||^. Figure 5 a shows the 
potential along the CCB(a)-UP direction (in the case Aq = 0, niQ = 100, mi/2 = 200, 
tan P{Mz) = 2, and /i < 0) using only the tree level potential. This point would have been 
excluded by the CCB(a)-UP constraint since V^hf' < Knin'^^- However, including the 1- 
loop correction in this case leads to the potential displayed in Fig. 5b. It appears unbounded 
from below as last seen at field values nearing the Planck scale in the tu-case thus also ruling 
out this point. However, in the same figure, the a-case 1-loop effective potential is not 
unbounded from below and indeed does not rule out this point (in this direction). 

IV. RESULTS 

Using the procedures outlined in Sec. HI and Appendices A and B, we explored regions 
of minimal SUGRA parameter space for minima deeper than the standard MSSM one. Our 
initial scans took place in the mo vs. mi/2 plane, to facilitate comparison with recent results 
on fine-tuning, cosmology and present and future collider searches. We fix tan f3{Mz) to 
be 2 or 10, take Aq = and rrit = 170 GeV. Our search was performed in the ranges 
< mo, mi/2 < 500 GeV and the grid was scanned with 25 GeV resolution. Figures 6a-d 
display the regions where non-standard global minima were discovered. Of all the directions 
scanned for these plots, non-standard vacua were found only in the UFB-3a direction. In 
Fig. 6, we have encoded information about the magnitude of the VEV in the plotting symbol. 
Using T] = logj^plf /wmssm}, the squares represent 2 < 77 < 3, the crosses 3 < t] < 4, and 
the x's 4 < ?7 < 5. The most dangerous regions are those populated by squares: For these 
points the "distance" between the standard and non-standard minima is smallest, which 
would admit the largest rate for tunnelling between them. Performing this scan using the 
exact prescription of Ref. leads to nearly identical excluded regions. 

We see from Fig. 6 that for all four frames, the region of low mo becomes excluded. 
As noted in Ref. [Q, this rules out the so-called "no-scale" models which require mo = 0. 
In addition, in string models where supersymmetry is broken in the dilaton sector, one is 
led to GUT or string scale soft-terms related by mi/2 = —Aq = v^mo [^,^. For this 
precise choice of soft-term boundary conditions, much of the parameter space is excluded 
by non-standard minima. We further note that the excluded region rules out much of the 
SUGRA parameter space associated with light sleptons. In particular, taken literally, our 
results exclude regions where such decays as Z2 —>■ ii^ and Z2 — *■ I'l' take place. 

Although Fig. 6 is plotted for Aq = 0, a similar excluded region results for other choices of 
the Aq parameter. This is shown in Fig. 7, where we plot regions excluded in the mo vs. Aq 
plane, for the same values of tan (3 and fi, but for mi/2 fixed at 200 GeV. The vacuum 
constraints exist for all Aq values, but are smallest for Aq ~ 300 GeV. 
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In Fig. 8, we display a combined plot of Fig. 6 with superposed dark matter and fine 
tuning [|ll[] contours. The regions to the right of the solid line contours are cosmologically 
excluded because they predict a relic density Qh'^ > 1; this corresponds to a lifetime of the 
universe of less than 10 billion years. Cosmological models which take into account COBE 
data, nucleosynthesis, and large-scale structure formation prefer an inflationary cosmology, 
with a matter content of the universe comprising 60% cold dark matter {e.g., neutralinos) , 
30% hot dark matter {e.g., neutrinos), and 10% baryonic matter. In this case, the preferred 
relic density of neutralinos should be .15 < flh'^ < .4, i.e., the region between the dot-dashed 
contours. In addition. Fig. 8 contains two naturalness contours with varying degrees of 



acceptability: 72 = 5 and 10 [|11[]. The more encompassing contour is a conservative estimate 
of a reasonable "tolerance limit" for weak scale supersymmetry. We see that the constraint 
from false vacua overlaps considerably with the preferred regions from cosmology and fine- 
tuning, leaving only a small preferred region of parameter space around mg ~ 100 — 200, and 
mi/2 ~ 100 — 250 in each frame. We note that the resulting preferred region of parameter 
space requires mg,mq,m^^,mi^ < 650, 600, 220 and 175 GeV, respectively, for tan/5 = 2, 
and mg,mq,m^^,mi^ ^ 725, 650, 225 and 220 GeV for tan/5 = 10. The only exception 
to these bounds is if the neutralino is poised near the peak of an s-channel pole in its 
annihilation cross section. These regions correspond to the narrow horizontal corridors in 
the relic density contours. 

Recently, the reach of the Fermilab Main Injector and TeV33 have been calculated in the 
same parameter space frames |13]. By comparing the results of Ref. [13] with the preferred 
parameter space discussed above, we see that the TeV33 option covers most of the preferred 
region from Fig. 8 a via the clean trilepton signal from W1Z2 3i, the exception being 
the region with mi/2 ~ 180 GeV, where the spoiler decay mode Z2 — > Zih turns on. For 
Fig. 86, TeV33 will cover the entire preferred region via clean trileptons. For the large 
tan/5 = 10 cases of Figs. 8c and d, the TeV33 upgrade can see most, but not all, of the 
preferred parameter space regions. The reach of the Tevatron Main Injector is significantly 
less than TeV33 for these preferred regions of parameter space. The CERN LHC collider 
can of course probe all the preferred regions of parameter space. In fact, event rates will 
be enormous for various multi-lepton + multi-jet +l^rp channels, which should facilitate 
precision measurements of parameters |T^. In particular, sleptons have mass less than 250 
GeV in these regions, and so ought to be visible at the LHC. Finally, we note that both the 
light chargino and right-selectron have mass less than 250 GeV in the preferred regions, so 
that both of these sparticles would be accessible to Next Linear Collider (NLC) experiments 
operating at i/i = 500 GeV [p 



V. CONCLUSION 

The minimal SUGRA model provides a well-motivated and phenomenologically viable 
picture of how weak scale supersymmetry might occur. The 4+1 dimensional parameter 
space can be constrained in numerous ways as discussed in the introduction. To constrain 
the model further, we have pursued the idea that parameter values that lead to global 
minima in non-standard directions, such as those with charge or color breaking, should be 
excluded from consideration. There are cosmological issues pertaining to tunnelling that 
we have knowingly ignored. Nevertheless, with this limitation noted, we searched for the 
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preferred regions of parameter space. We analyzed the potentials carefully employing the 
1-loop correction including the contributions of the top and stops in the calculation. In 
generating the potential we found that the 1-loop correction can significantly alter results 
based only on the tree approximation. 

Because of the various scales present in the MSSM, renormalization group improvement 
has ambiguities associated with it. We tried several procedures but were ultimately led to 
two distinct results that we refer to as the a- and u;-cases. In the cu-case, to our surprise, 
the entire parameter space of the model suffers from global minima along non-standard 
directions (namely the UFB-3(b) and CCB(a)-UP directions). Since the 2-loop correction 
may be significant, given the large value of the top mass, the results in the a;-case, while 
intriguing, must be taken cum grano salis. 

In the a-case, we are still left with a very restricted region of parameter space after 
imposing in addition dark matter and naturalness constraints. Most of this region should 
be accessible to the Fermilab TeV33 collider upgrade via the clean trilepton channel. This 
parameter space region should be entirely explorable at the LHC, and should yield a rich 
harvest of multilepton signals for supersymmetry which ought to allow for precision deter- 
mination of underlying parameters. In addition, both charginos and sleptons ought to be 
accessible to NLC experiments operating at just ^/s = 500 GeV, so that the underlying 
assumptions of the minimal SUGRA model can be well-tested. 
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APPENDIX A: 



In this appendix, the cases from that are considered in our investigation are reviewed. 
Also, formulas for the top/stops are displayed in each case and for the bottom/sbottoms in 
the last case. 

• UFB-1: In the UFB-1 case, the only fields acquiring non-zero VEVs are < >=< 
$d >= X. The resulting tree level potential in this direction is 



V = (jnl + m\ — 2|m3|^ 



(A.l) 



where m\ = + /x^, m\ = + /i^, and 
stop mass matrix entries are 



rrin 



Bfi. The top mass is Mt = ytx, and the 



M 



LL 



2,22 



yt {fi - At)x . 



(A.2) 



• UFB-2: In the UFB-2 case, an additional slepton field is included to help control the 
D-terms. The shifted fields are therefore 



<^u> 

< $d > 
< L3 >^ 



X , 



(A.3) 



The z/ subscript represents the SU(2) direction that has acquired the VEV. The scalar 
potential in this case is 



2 22, 22 r,|2|2| 2 2 2 . ^ "2 

V = m^x + rrigX — 2,'-)\m.^\x + 7^m^^x + ^fl' 



7 -7l 



X 



(A.4) 



Minimization with respect to 7 and 7I gives 



7 



2 2 ' 
mf — m|^g 



f - 



2mi3 
g'^x'^ 



(A.5) 
(A.6) 



If 7I < 0, then 7^ = 0, and we recover the UFB-1 direction. In this case, the top mass is 
Mt = ytX, and the stop mass matrix entries are 



^Q, + Vt X + [ T7;9 - 7^2 1 - 7 - 7l 



1 



,'2 



12 



2,22 ^ 

+ yt X - 3^ 



/2 



4^ 



1 2 2 

1-7 -7l 



X 



X 



^LR - yt (/^7 -A)x. 

UFB-3(a): In the UFB-3(a) case, the shifted fields are 



(A.7) 

(A.8) 
(A.9) 



12 



< > = a; , 



.fjLX . 



Vr 

< L2>u = Ilx . 



(A.IO) 



The VEVs of the (63, 63) sleptons are fixed by imposing $d F-term cancellation in the 
potential. The scalar potential for this case is 



rrir- 



Minimization with respect to 7^ gives 



1 r i2 
2, 2 2 2|-'--2 2i/)2 22 



(A.ll) 



72 = 1 + 



2mL 



(A.12) 



If 7^ < 0, then ^\ — 0. In this case the top mass is Mf — i/tx, and the stop mass matrix 
entries are 



1 

12' 



J2 



X -\- 1 — 72,x 



Ti/r2 2,2 2 ^ II 1 . 1)2 22 



MIr = -Atvtx . 
UFB-3(b): In the UFB-3(b) case, the shifted fields are 

<^d> = , 

< Q3 >d = < 4 = 

< L3 >^ = 'JlX . 



(A.13) 

(A.14) 
(A.15) 



fj,x I 

Z/6 



(A. 16) 



The VEVs of the (^3, ^3) squarks are fixed by imposing F-term cancellation in the 
potential. The scalar potential appears as 



jr 2 2 , ( 2 I 2\j2i 2 2 2|-'--2r2|j2 2' 

V = m^^x + [niQ^ + m^J + 7^771^32; + -c/ + d - 7^2; 
Minimization with respect to 7^ gives 



(A. 17) 



72 = 1 + 



2m|3 
^2^2 



(A.18) 



If 7^ < 0, then 7^ = 0. In this case the top mass is Mt — ytx, and the stop mass matrix 
entries are 



13 



M-2 2 I 2 2 I 2 j2 I '2 r 2 I j2 2 2 

LL = ^Q, +yt X +yb d + T7;9 X +d --f^x 



1 

12' 
1 



1 2 
4^2 



2 72 2 2 

X — a — 



71 /r2 2 I 2 2 I 2 j2 /2 2 , j2 2 2 



CCB(a)-UP: In the CCB(a) case for the up-trihnear, the shifted fields are 

<^d> = , 
<Qi>u = ax , 
<Ui> — Px , 

<Qs >l ^<ds>'^d'^\^\ , 



Vb 



(A.19) 

(A.20) 
(A.21) 



(A.22) 



SU(3) D-flatncss impHcs /?2 = a^. Also, U(l) and SU(2) D-flatness imply 1 - - jI = 0. 
Consequently, the scalar potential appears as 

V^yl(2 + q;2) a'^x^ - 2Tia^x^ + (m'^o^ + + mfj x^ , (A.23) 



where Ti = |A„|/„a;|, 
q;^ yields 



Qs + "^is ~ = Minimization with respect to 



, _ Ti - MV2 - ygx^ 

9 9 



(A.24) 



If q;^ < 0, then = and 7I = 1. If 7I < 0, then one should try < L3 >e=< 63 >= JlX 
and < L3 >u= 0. D-flatness now implies 1 — a2 + 7|, = 0. This also changes mj = m\^+m\^. 
In this case the top mass is Mt = ytx, and the stop mass matrix entries are 



Ml, 



m' + y1 x^ + yi + ( ^-g''' + -a^ ) d' , 



.2 j2 



'2 



12^ 



M 



B.R 



2 I 2 2 I 2 j2 /2j2 

^ua + Vt X +yt d - -g d , 



CCB(b)-UP: In the CCB(b) case for the up-trilinear, the shifted fields are 

< > = X , 

< $d > = 7a: , 

< Qi >« = ax , 

< Ml > = Px , 

< L3 >^ = 7La; . 



(A.25) 

(A.26) 
(A.27) 



(A.28) 



Unlike the previous CCB case (7 = 0), this case has three terms whose phases (q = cos (fi) 
are undetermined 
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2\AuyuUi'^uQi\ci + 2\fiyuUi<5lQi\c2 + 2|/i5$„$rf|c3 



(A.29) 



Reference |^ shows that this ambiguity can be resolved into two distinct possibihties. If 
sign(y4u) = — sign(i?), the three terms can be made simultaneously negative. If sign(y4„) = 
sign(i?), then the term of smallest magnitude is taken positive and the other two can be 
taken negative. 



SU(3) D-flatness implies = a^. Also, U(l) and SU(2) D-flatness imply 1 



72 = 0. Consequently, the scalar potential appears as 

V = yl{2 + a^) aV + 2 (ciTia^ ^ c2T2a^\i\ + C3r3|7 



X 



+ 



M^a^ + {ml - mj)Y + m^ + mj 



X 



where 











Ts 




M2 


2 1 2 

= mr, + m^ — 


ml 


2 1 2 


ml 


2 1 2 



^1 ! 



a2 - 72 - 



(A.30) 

(A.31) 
(A.32) 
(A.33) 
(A.34) 
(A.35) 
(A.36) 
(A.37) 



and where m| = m^^ and q = cosy^j. Minimization with respect to and I7I (note that 
the potential is a function of I7I) yields 

02X2^2 + C3T3 



171 



a 



2 2 



i/gx2 + MV2 + ciri + c2r2| 



(A.38) 
(A.39) 



It must be confirmed that both > Q and I7I > 0. Otherwise these are set to zero. It 
must also be checked that 7^ > 0. Otherwise one should try < L3 >e=< 63 >= 7lx and 
< L3 >u= 0. D-flatness now implies 1 — a2— ^2_j_^2 _ g_ ^jj-^jg ^[^q changes ml = m\^+m^ 
In this case the top mass is Mj = ytx, and the stop mass matrix entries are 



63' 



Ml, 



M: 



RR 
2 

LR 



2,22 

^Q, + yt X , 

2,22 
^u, + yt X , 

yt (/i7 - At)x . 



(A.40) 
(A.41) 
(A.42) 



• CCB(b)-TOP: In the CCB(b) case for the top-trilinear (there is no (a) case), the shifted 
fields are 



< $a > 

< $d > 

< Q3 >u 
<U3> 

< L3 >^ 



X , 

7X , 
ax , 
(3x , 



(A.43) 



15 



SU(3) D-flatncss implies 0"^ = o?. There is no imposition of D-flatness in the U(l) and 
SU(2) sectors. The potential now appears as 



2/01 2\24ir2/2 I 2\| 2 2, 2 2 , 2 

^ = ?/t (2 + a )a X + [a {mq^ + m^g) + 7 + 7^771^3 + 



X 



+ 2 (ciTia^ + c2T20?Yi\ + C3T3I7I) + hj" (l 



2 2 2\ 4 

a -7 -7Lj^ 



(A.44) 



where 





= \Aytx\ , 


T2 


= \yt^J'X\ , 


T3 






= 772^3 + - m% 




2 1 2 




2 1 2 



(A.45) 
(A.46) 
(A.47) 
(A.48) 
(A.49) 
(A.50) 



Minimization with respect to 7^ yields 



2 1 22 2771? 
7l = 1 - a - 7 - 

g-^x^ 



If 7^ > 0, then minimization with respect to and I7I gives 

I I C2T2a^ + C3T3 



m 
,2^2 



2 



|/,V + MV2 + ciTi + C2T2|7| 



(A.51) 

(A.52) 
(A.53) 



If 7I < 0, then it must be set to zero, and minimization with respect to and I7I yields 

2y,V (1 + a^) - ^g^x' (1 - - 7') + {ml„ + n^^) x^ + 2 (ciTi + C2T2I7I) = (A.54) 
-^V (1 - a' - 7^) I7I + 2|7|xV + 2 (c2r2a2 + C3r3) x^ ^ . (A.55) 

Substituting for I7I yields a cubic equation for a^. It must still be checked that both I7I > 
and > 0. In this case the top mass is Mt — i/tx, and the stop mass matrix entries are 



ML 



LR 



ml.^+ylx\l + a') + —g'^ 



1^2 2 2 
1-3" - 7 - 7l 



X 



:9l 



1 



1 o2 2 2 2,-^2 

1 - 3q; - 7 - 7^J a; + -^3 



2 2 
a X 



2 



2 I 2 2/1 I 2\ /2 



1 ' 2 2 2 

1-3" - 7 - 7l 



2 I 2 
2; +3^3 



2 2 



l/i (/i7 - At) X + 



1 /I 



- 3 ^3^ 



x^a(5 . 



(A.56) 
(A.57) 
(A.58) 



• CCB(a)-ELECTRON: In the CCB(a) case for the electron-trihnear, the shifted fields 
are 



16 



< Li >e = ax , 
< ei > — Px , 



<Q3>'t ^<U3 U' = 



I A*^ I 
Vt 



(A.59) 



D-flatness implies = and a^x^ = + m^. The scalar potential appears as 
V^yl{2 + a^) a'^x^ + [m|^ + Q;^(mi^ + m^j] + (m^g + m^^) - 2Tro?x^ , (A.60) 



where Ti = |Ae|/ea^|- In this case we use the bottom/sbottom contribution. The bottom 
mass is = yhX, and the sbottom mass matrix entries are 



MIl = ml + yy 



1 



Vt 



;92 u 



I + yl + ' 



(A.61) 

(A.62) 
(A.63) 



• CCB(b)-ELECTRON: In the CCB(b) case for the electron-trihnear, the shifted fields 
are 



< *d > 

< > 

< ei > 



x 



ax , 
j3x . 



(A.64) 



D-flatness again implies 0^ — o? and = 1 — 7^. The potential appears as 



where 



yl (2 + a") a'x' + 




+ 2 (ciTia^ + C2T2 


«^|7l + C3T3I7I) a;^ , 


Ti = 1 




T2 = | 


y^fix 1 , 


T3 = \^^B\ . 



2^ 



(A.65) 



(A.66) 
(A.67) 
(A.68) 

Substituting for and minimizing with respect to I7I gives the following cubic equation 



|7| 



2yy - f + |7| (ml-ml^-ml)-4yy-2ciT, + [caTs + C3T3] - 



(A.69) 



It must be checked that both I7I > and a^ > 0. In this case the top mass is Mt — "yytx, 
and the stop mass matrix entries are 
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= rn%^ + lYt , (A.70) 

MIr = + -iVt , (A.71) 

Mlj,^yt{lx-^At)x . (A.72) 

If 7 = 0, then we use the bottom/sbottom contribution. The bottom mass is = y^a;, and 
the sbottom mass matrix entries are 

MIl = rnl^ + Vl x" , (A.73) 

MIr = A^yftX . (A.75) 
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APPENDIX B: 



We find that Eq. (2.3) cannot be reliably calculated using our computer for large (> 
10^) values of the VEVs. We therefore use a limiting form of this expression in the large 
VEV hmit. This hmiting form of AVi is presented in this appendix. We begin with some 
definitions: 



mt = ytx 

= W + + Mlx + Glx^ 
^RR = + + ^R^ + Grx^ 



(B.l) 
(B.2) 
(B.3) 
(B.4) 



These expressions cover all the cases we have analyzed, and it is a simple matter to identify 
the various coefficients for each case. Note that in the CCB-ELECTRON cases, rrit should 
be substituted with m^. In terms of these definitions, the two stop masses are 



t± 2 
± 



2m^ + {ml + m\) + {Ml + Mr)x + {Gl + Gr)x^ 



[ml - ml) + {Ml - Mr) + {Gl - Gr)x' 



+ 4 



Mx + Gx' 



1/2 



(B.5) 



To simpifiy the notation some new definitions are used, and the stop masses rewritten in 
terms of these 



where 



Ml = m', 



± -m. 




mt 



r + —r + 

mi 



—0 + — e 



1/2 



(B.6) 



G± 
M± 
m\ 
a" 

,2 



Gl ± Gr 

Ml±Mr 

2 I 2 

± 

{G^ + 4G')/yt 
2(G_M_ + 4MG)/y^ 
{Ml + 2mlG- + 4:M^)/yt 
2mlM_/yt 

4 

m . 



(B.7) 
(B.8) 
(B.9) 
(B.IO) 
(B.ll) 
[B.12) 
[B.13) 
(B.14) 



Modulo overall factors the 1-loop correction is 



AV^i = -2m^ln{m,7Q'} + Mf ln{M/ /Q'} + ln{M^/g^} 



(B.15) 



where Q = Qe'l'^ . There are three cases to be considered 
• case (a): o? ^ 0, 
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case (b): = q (=> /?2 = 0); 7^ 7^ 0, 





1 + 




B± = 


M+ 


±^ 
2 






±? 

47 



case (c): ^ ^ ^ ^ 0; ^ 0, 









1 + 






'2yt 




M+ 














2 






±'- 




2 



Finally the form of the 1-loop correction in the large rrit limit is 

AT4 = [2(a+ + a_) + (a^ + a^)] L + (^^ + A'^_ In A_) 

+ — U+B+(l/2 + L + In A+) + A_5_(l/2 + L + In 
mt 

+ ^[5^(3/2 + L + ln^+) + 5!(3/2 + L + ln^_) 

+ 2^+C+(l/2 + L + ln/l+) + 2A_C_(l/2 + L + \nA_)]} (B.25) 

where L = ln{mj/(5^} and a± = A± — 1. To arrive at AVi, multiply the above expression 
by 3/327r2. 
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FIGURES 



FIG. 1. Plots of the 1-loop correction, tree and 1-loop effective potentials along the MSSM and 
UFB-3(a) vacuum directions in the = 0, mo = 100 GeV, mi/2 = 200 GeV, /x < 0, and tan /3 = 2 
case. Renormalization group improvement was implemented using the a- prescription. 

FIG. 2. Plots similar to Fig. 1 but with renormalization group improvement implemented using 
the ^-prescription. 

FIG. 3. Plots of the logarithm of the VEV versus x-, ^ parameter appearing in the renormal- 
ization group improvement function, for the vacuum directions UFB-3(b) and CCB(a)-UP. 

FIG. 4. Evolution of the VEV by using the renormalization group 7 functions of the Higgs 
fields (solid line) and by tracking the minimum of 1-loop potential (dashed line) as a function of 
the renormalization scale Q. 

FIG. 5. Plots of the potential in the CCB(a)-UP case. In (a) the tree potential is displayed. In 
(b) the 1-loop potential is displayed for both a and uo precriptions. 

FIG. 6. Exclusion plots for the tuq vs. mi/2 plane based on the UFB-3(a) constraint. The 
squares represent 2 < 77 < 3, the crosses 3 < 77 < 4, and the x's 4 < 77 < 5. 

FIG. 7. Exclusion plots for the mo vs. Ao plane based on the UFB-3(a) constraint and with 
mi/2 = 200 GeV. The symbols are the same as in Fig. 6. 

FIG. 8. Same as Fig. 6 but with superposed dark matter (dot dashes and solid) and naturalness 
(dashes) contours. 
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